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We describe measurements of spin-polarized tunnelling via discrete energy levels of single Alu- 
minum grains. In high resistance samples (~ GO,), the spin-polarized tunnelling current rapidly 
saturates as a function of the bias voltage. This indicates that spin-polarized current is carried only 
via the ground state and the few lowest in energy excited states of the grain. At the saturation 
voltage, the spin-relaxation rate T-j -1 of the highest excited states is comparable to the electron tun- 
nelling rate: T^ 1 ~ 1.5 • 10 6 s~ 1 and 10 7 s _1 in two samples. The ratio of to the electron-phonon 
relaxation rate is in agreement with the Elliot- Yafet scaling, an evidence that spin-relaxation in Al 
grains is governed by the spin-orbit interaction. 



INTRODUCTION 

Electron tunnelling through single nanometer scale 
metallic grains at low temperatures can display a dis- 
crete energy level spectrum, [lj Tunnelling spectroscopy 
of the energy spectra have led to numerous discover- 
ies, including Fermi-Liquid coupling constants between 
quasiparticles, Q spin-orbit interactions, [1,0] and super- 
conducting correlations in zero-dimensional systems. [BJ 
Some information regarding the spin of an electron oc- 
cupying a discrete level can be obtained using spin- 
unpolarizcd tunnelling, such as spin-multiplicity and 
electron g-factors. [l| 

In this letter we report on spin-polarized tunnelling via 
discrete energy levels of single aluminum grains. Spin- 
polarized electron transport permits studies of spin re- 
laxation and spin dephasing. 0] By comparison, spin- 
unpolarized spectroscopy is suitable for the studies of en- 
ergy relaxation in the grains. [l|, Q Since spin-relaxation 
times are generally many orders of magnitude longer than 
energy relaxation times, spin-unpolarized spectroscopy is 
not an easy tool to study spin-relaxation in the grains and 
spin-polarized tunnelling is needed. We find that some 
electron spin-relaxation times in Al grains are exception- 
ally long compared to bulk Al with similar disorder, on 
the order of /is. 

Spin-polarized transport via metallic grains has re- 
cently generated a lot of theoretical interest. 



In addition, there is a major effort to study nano- 
spintronics using carbon-nanotubes; see Ref. [l3T ] and ref- 
erences therein. Spin-coherent electron tunnelling via 
nanometer scale normal metallic grains has been con- 
firmed in arrays [HE! and in single grains. 3] How- 
ever, the electron spin-relaxation time T\ in a metallic 
grain has not been reported yet. 



SAMPLE FABRICATION 

Our samples are prepared by electron beam lithogra- 
phy and shadow evaporation, similar to the technique 
described previously. [1| First we define a resist bridge 
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FIG. 1: A, B, C: Sample fabrication steps. D. Image of Al 
grains. E, F: Image of a typical sample. G, H: I-V curves at 
the base temperature. 



placed 250 nm above the Si wafer; this bridge acts as a 
mask. Next (Fig.[T]-A), we deposit 11 nm permalloy (Py 
= Nio.sFeo.2) onto oxidized silicon substrate at 4 • 10~ 7 
Torr base pressure, measured near the gate valve, along 
the direction indicated by the arrow. Then we rotate the 
sample by 36 degrees without breaking the vacuum and 
deposit 1.2 nm of AI2O3 by reactive evaporation of Al, Q 
at a rate of 0.35 nm/s, at an oxygen pressure of 2.5 • 10~ 5 
Torr. Now, oxygen flow is shut down. When pressure 
decreases to the 10 Torr range, we deposit a 0.6 nm 
thick film of Al, as shown in Fig. [l]-B. Al forms isolated 
grains with a typical diameter of 5 nm. The grains are 
displayed by the scanning electron microscope (SEM) im- 
age in Fig. [l]-D. Finally we deposit another 1.2 nm layer 
of AI2O3 by the reactive evaporation and top it of by 
an 11 nm thick film of Py (Fig. [T]-C). We make many 
samples on the same silicon wafer, and vary the overlap 
from to 50 nm and select the devices with the highest 
resistance, as they have the smallest overlap. Figs. [T]-E 
and F show SEM images of a typical device. 
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DISCRETE ENERGY LEVELS 

Transport properties of the samples at low tempera- 
tures were measured using an Ithaco current amplifier. 
The samples were cooled down to « 0.035K base tem- 
perature. The sample leads were cryogenically filtered to 
reduce the electron temperature down to f=a 0.1K. 

The majority of samples (>80%) exhibit Coulomb 
Blockade at low temperature. About 150 samples were 
measured at 4.2K and 16 samples at 0.035K . In this pa- 
per we describe two samples. The I-V curve of two sam- 
ples are shown in Fig.[l]-G and H. The tunnelling current 
increases in discrete steps as a function of bias voltage, 
corresponding to discrete electron-in-a-box energy levels 
of the grain. 

In sample 1, the average electron-in-a-box level spacing 
caused by electron geometric confinement is 5 ~ 0.8meV, 
which corresponds to diameter of D w 6nm assuming 
a spherical Al grain. The average current step is I w 
0A7pA. We make a connection with the tunnelling rates 
from the leads to the grain and the measured current 
response. The tunnel junctions are highly asymmetric, 
and therefore one of the tunnelling rates is much smaller 
than the other, and thus rate limiting. Throughout this 
paper, we choose the rate limiting step to be across the 
left junction, corresponding to the tunnelling rate T^. 
Therefore, our measured current corresponds to the av- 
erage tunnelling-in rate of Tl = Z/2|e| « 1.5 ■ 10 6 s _1 . 
Similarly, in sample 2, 5 k 2.7r7ieV, D w 4nm, and 

r L wg.e-io 8 *- 1 . 

The spin-conserving energy-relaxation in Al grains 
takesplace by phonon emission with the relaxation 
rate [2j 
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where Ep = ll.7eV is the Fermi energy, uj is the energy 
difference between the initial and the final state, p = 
2.7 'g I 'cm 3 is the ion- mass density, and v s — 6420m/s is 
the sound velocity. We obtain T^} ph (5) w 1.6 • 10 9 s _1 
and 4.1 • 10 10 s _1 in samples 1 and 2, respectively. Sample 
2 has significantly larger relaxation rate because of the 
larger level spacing. Since the tunnelling rates in our 
samples are ~ 10 6 s _1 , if the grain is excited by electron 
tunnelling in and out, it will instantly relax to the lowest 
energy state accessible by spin-conserving transitions. 

As shown by Fig. [21 the energy levels exhibit Zeeman 
splitting as a function of an applied magnetic field. In 
sample 1, the I-V curve probes the same energy spec- 
trum at negative and positive bias voltage. This is evi- 
dent from the equivalence of the magnetic field dependen- 
cies at negative and positive bias. The lowest tunnelling 
threshold is two fold degenerate at zero magnetic field, 
showing that Nq, the number of electrons on the grain 
before tunnelling in, is even. The conductance peaks are 
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FIG. 2: A, B: Differential conductrance (gray) versus bias 
voltage and the applied magnetic field in sample 1 at the 
base temperature. 



similar in magnitude at negative bias, because the first 
tunnelling step, in which an electron tunnels in to the 
grain through the higher resistance junction, is rate lim- 
iting. At positive bias, the first conductance peak is much 
larger than the subsequent conductance peaks, because 
the first tunnelling step takes place via the lower resis- 
tance junction, and the rates are limited by the electron 
discharge process across the high resistance junction. 

In sample 1, the first two peaks split corresponding to 
g-factors: g = 1.83±0.05 and 1.95±0.05. Slight reduction 
of the g-factors from 2 indicates spin-orbit interaction in 
Al. [l| The avoided level crossings clearly are resolved in 
Fig.[l near points (-11.5mV,5T) and (-13mV, 11.5T). 
The corresponding avoided crossings at positive bias are 
located near (13.5mV, 5T) and (15.5ml/, 11. 5T), respec- 
tively. In the regime where g factors are slightly reduced, 
the spin-orbit scattering rate (t^q) can be obtained from 
the avoided crossing energies Ago ~ 0-lmeV. [17| The- 
ory predicts that t$o ~ M/ttAjq, 
tor of two. Thus, we obtain t^q ~ 

" T SO 5 

= OLT. 



17] within a fac- 



5.5 ■ 10 w s~ l . By 
is related to the elas- 
Assuming ballis- 



10 14 sec 1 . We obtain 
in excellent agreement with a ~ 10~ 4 in 



the Elliot- Yafet relation, [18( 
tic scattering rate T~ l : t^q 
tic grain, r" 1 « vp/D = 3.4 
a w 1.6 • 10"' 
Al thin films. 



SPIN-POLARIZED TUNNELLING 

Now we discuss magnetoresistance from the spin- 
polarized tunnelling. In the magnetic field range of 
±50mT, approximately 90% of the samples do not 
display any of the tunnelling magnetoresistance effect 
(TMR) . By contrast, we tested about 10 tunnelling junc- 
tions without the embedded grains and with similar re- 
sistance (empty junctions) at 4.2 K. All of the empty 
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junctions exhibit a significant TMR in this field range, 
comparable to 10%. Approximately one half of the empty 
junctions display a simple spin-valve effect. So, the ab- 
sence of TMR for electron tunnelling via grains shows 
that the spin-dephasing rate T^ 1 in 90% of the samples 
must be much larger than the tunnelling rate. 

Nevertheless, approximately 10% of the samples with 
embedded grains display significant TMR, so the dephas- 
ing must be weak, e.g. T 2 _1 must be smaller than or com- 
parable to the tunnelling rate in these samples. Ti varia- 
tion among different samples could be explained by mag- 
netic defects, such as paramagnetic impurities from the 
Py layer. Paramagnetic impurities are common sources 
of dephasing. The defects would be located on the 
grain surface, since bulk Al does not support paramag- 
netism. Since the number of atoms on the surface is 
relatively small (~1000), we could occasionally obtain a 
sample free of impurities. More insight into the nature of 
T% in this device will require a more in depth theoretical 
study. 

A majority of the samples with nonzero TMR show 
positive TMR near the Coulomb-Blockade conduction 
threshold; only about 30% of the samples show nega- 
tive TMR. The sign of TMR in quantum dots is deter- 
mined by the inter play between charging effects and spin- 
accumulation. TH, 21] For any given sample, the data in 
this paper correspond to the voltage range within the first 
step of the Coulomb staircase. In this range the sign of 
TMR is found to be constant as expected. 

TMR in our devices usually does not display a simple 
spin-valve effect. We believe this is because there are 
spin-dependent interactions inside the grain that induce 
a complicated TMR even when the magnetic transitions 
in the drain and source leads are sharp and as expected. 
For example, a rotation of stray magnetic field acting on 
the grain will alter the direction of the spin-quantization 
axis in the grain, thereby changing the conductance [8j]. 
A rotation or a switch of a remote domain can change 
the tunnelling current through the grain via the mag- 
netic field generated by the domain. Similarly, the ori- 
entation of the nuclear spin in the grain can change the 
quantization axes via the hyperfine interaction. 

We select only those samples that display a simple 
spin- valve TMR effect, which is shown in Figs. [3] Fig. [3J- 
A is the TMR of sample 1 at a bias voltage corresponding 
to the second current plateau. TMR is barely resolved in 
this case, since the current changes by only about 40fA. 
We do not have good data to display TMR at the first 
current plateau. By comparison, Figs. [3]- B and C display 
TMR at bias voltage where the number of electron-in- 
a-box levels energetically available for tunnelling-in are 
approximately 19 and 48, respectively. To facilitate com- 
parisons, the current intervals on the vertical axes in 
Figs.[3]A-C and D-F have equal lengths. 

The main observation in this letter is that AI = 
7jl —Ifl is nearly constant with current above a cer- 
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FIG. 3: A-F: Spin- valve effect in current versus applied mag- 
netic field in two samples at the base temperature. The cur- 
rent magnitude is reduced in the antiparallel state. 



tain current. There is hardly any increase in AI between 
Figs. [3] B and C and between Figs. [3] E and F. This be- 
havior is shown in more detail Fig. |4]A and B, which 
displays AI versus bias voltage. AI versus negative bias 
voltage in sample 1 is fully saturated at the third current 
plateau; at the second current plateau, AI is already at 
one half of the saturation value. Similarly, in sample 
2 AI reaches saturation at the second current plateau. 
Our samples should be contrasted with ordinary ferro- 
magnetic tunnelling junctions, where AI is proportional 
to the current over a significantly wider range of bias 

voltage (mm. 



INTERPRETATION OF THE RESULTS 

In coulomb-blockade samples containing magnetic 
leads, the electrochemical potential difference between 
the island and leads can jump when the magnetization 
in one of the leads changes direction [1(3]. This can lead 
to a sudden shift in energy levels, producing a jump in 
current that is constant as a function of bias voltage. The 
shift in energy levels is seen as a discontinuity near zero 
magnetic field in Fig. 2, and is ~ 0.1rnV\ 

To show that the electrocemical shift is not responsi- 
ble for the saturation of the spin-polarized current with 
voltage in our sample, we performed other measurements 
by sweeping the magnetic field both on and between the 
current plateaus, coming up with similar values for the 
electrochemical shift. The shift is lower than the aver- 
age level spacing of 0.8 meV and 2.7 meV for sample 1 
and sample 2 respectively. Therefore, since we measured 
magnetoresistance in the middle of the current plateau, 
the threshold voltage shift should not effect our measure- 
ments of the saturation in AI. 

To explain iff — Iff = const, we must discuss the rel- 
ative magnitudes of three rates: r~_ h , the rate of en- 



ergy relaxation from excited to lower energy states by 
spin-conserving phonon emission; T^, the rate electrons 
tunnel into the grain; and T^ 1 , the rate of transitions 
between levels that result in an electron flipping its spin 
orientation. T~_ ph is obtained theoretically, the measured 
I-V spectrum fixes the tunnelling rate, and Tj -1 is ob- 
tained from the saturation in J-j-j- — J-j-j with bias voltage. 

Finally we must deduce the relative magnitude of I\ . 
The rate of spin-flip transitions is expected to be signif- 
icantly smaller than t~_ h |18j. In this case the ground 
state would not necessarily be accessible by energy re- 
laxation. The grain could remain in an excited, spin- 
polarized state, as sketched in Fig. 0J-C. These spin- 
polarized excited states are responsible for spin accumu- 
lation in the antiparallel magnetic configuration of the 
leads. If the relaxation rates for the spin-flip transi- 
tions arc much smaller than the tunnelling rate, then 
various spin-polarized states would have similar proba- 
bilities, which are determined by the tunnelling rates. In 
the antiparallel configuration of the leads, the probabil- 
ities of the excitations with spin up would be enhanced 
by 1 + P and probabilities of the excitations with spin 
down would be suppressed by 1 — P, where P is the spin- 
polarization in the leads. In the parallel configurations, 
the probabilities of the excitations with spin up and spin 
down are the same. In this regime, iff — iff is propor- 
tional to the current, similar to the usual ferromagnetic 
tunnelling junctions. 

It is reasonable to expect that the spin-flip rate T 1 _1 (w) 
increases rapidly with e nerg y difference uj between the 
initial and the final state 24]. If T 1 _1 (w) exceeds the tun- 
nelling rate above some uj, then the excitations with en- 
ergy > u) will occur with a reduced probability in the en- 
semble of states generated by tunnelling in and out. Thus 
AI is limited by tunnelling via the ground state and those 
low lying spin-polarized states where T 1 _1 (w) < T^. AI 
versus bias voltage approaches saturation approximately 
when T^ 1 {uj) — Tl, where uj is the highest excitation 
energy in the ensemble of spin-polarized states generated 

This is how we 



Telr7 



by tunnelling in and out: uj 
determine the spin-relaxation time T\ (uj) at an energy uj 
in a given sample. 

In sample 1, AI is at 50% of the saturation value at the 
second current plateau, and AI is saturated at the third 
current plateau. At the second current plateau, the spin- 
relaxation rate of the highest energy excited state gener- 
ated by tunnelling must be close to the tunnelling rate. 
Since the spin relaxation is very rapid in configurations 
more than 3(5 above the ground state, and Ao is even as 
noted above, the grain spends most of the time among the 
five configurations shown in Fig.0]-C: Nq, Nq, Nq , N^ + , 
and Nq~. The highest energy spin-polarized states are 
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Thus, Tf L (3S) w T L = 1.5 ■ 10 ( 
sample 2, this analysis leads to T^ 1 (2S) r* 10 7 s" 
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FIG. 4: A and B: AI = 



1 1l\ versus bias voltage in 



samples 1 and 2, respectively, at the base temperature. The 
numbers near the circles indicate how many doubly degen- 
erate electron- in-a-box levels are available for tunnelling in. 
C: Possible spin-polarized electron configurations caused by 
electron tunnelling in and out, before an electron tunnels in, 
at the second current plateau, for Nq even. 



Now we discuss the origin of spin relaxation and its 



rapid enhancement with the energy difference. Note that 
the rate of spin-conserving transitions in Eq. [1] increases 
as uj 3 . We suggest that the electron-phonon transition 
rates without and with spin flip scale by the Elliot- Yafet 
relation: T 1 _1 (cj) = a'T~_ ph (oj). This scaling would cer- 
tainly explain the rapid increase in spin-relaxation rate 
with excitation energy. In metallic films, it is well estab- 
lished that the Elliot- Yafet scaling applies for both elastic 
and inelastic scattering processes, with a rj a'. (l9j 



In sample 1, Eq. Q] leads to T~* ph (3S) » 4 • 10 10 s _1 . 
Since T^iSS) w 1.5 • 10 6 s-\ we obtain a' « 0.4 • 10~ 4 . 
Similarly, in sample 2, t~_ h (2S) ~ 3.3 • 10 11 s _1 and we 
obtain a' « 0.3 • 10~ 4 . a' agrees with a « 1.5 • 10~ 4 
obtained earlier, within an order of magnitude. So the 
ratio of T e - p h and T\ is in agreement with the Elliot- 
Yafet scaling. This is an evidence that the spin-flip tran- 
sitions in AI grains are driven by the spin-orbit interac- 
tion. By this relaxation mechanism, the spin of an elec- 
tron on the grain is coupled to the phonon continuum 
via the spin-orbit interaction. An electron in an excited 
spin-polarized state relaxes by an emission of a phonon, 
which has an angular momentum equal to the difference 
between the initial and final electron spin. 
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CONCLUSION 

In summary, we have observed spin-coherent electron 
transport via discrete energy levels of single Al grains. 
Spin polarized current saturates quickly as a function of 
bias voltage, which demonstrates that the ground state 
and the lowest excited states carry spin polarized current. 
Higher excited states have a relaxation time shorter than 
the tunnelling time and they do not carry spin-polarized 
current. The spin- relaxation time of the low- lying excited 
states is 7\ « 0.7 /is and 0.1/is in two samples. Finally, 
the ratio of the spin-flip transition rate and the electron- 
phonon relaxation rate is in quantitative agreement with 
the Elliot- Yafet scaling ratio, an evidence that the spin- 
relaxation transitions are driven by the spin-orbit inter- 
action. 
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